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A 2D electron system in a quantized magnetic field can be driven by microwave radiation into a 
non-equilibrium state with strong magnetooscillations of the dissipative conductivity. We demon- 
strate that in such system a negative conductivity can coexist with a positive diffusion coefficient. 
In a finite system, solution of coupled electrostatic and linear transport problems shows that the 
diffusion can stabilize a state with negative conductivity. Specifically, this happens when the system 
size is smaller than the absolute value of the non-equilibrium screening length that diverges at the 
point where the conductivity changes sign. We predict that a negative resistance can be measured 
in such a state. Further, for a non-zero difference between the work functions of two contacts, we 
explore the distribution of the electrostatic potential and of the electron density in the sample. 
We show that in the diffusion-stabilized regime of negative conductivity the system splits into two 
regions with opposite directions of electric field. This effect is a precursor of the domain structure 
that has been predicted to emerge spontaneously in the microwave-induced zero-resistance states. 
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I. INTRODUCTION 



High mobility two-dimensional electron systems 
(2DES) subjected to the microwave radiation reveal gi- 
ant magnetooscillations in the diagonal resistance Pxx^^ 
with periodicity controlled by the positions of the cy- 
clotron resonance harmonics. In high-quality samples, 
these oscillations give rise to zero resistance states^^^ 
(ZRS) at the oscillation minima, pxx -^ 0, when the ra- 
diation power is increased and the temperature is low- 
ered. In the Corbino-disk geometry ZRS manifest them- 
selves as zero-conductance states. The widely accepted 
explanation of ZRS combines two distinct phenomena. 
The first one is development of the giant magnetoresis- 
tance oscillations due to indirect optical transitions^"— or 
nonequiUbrium occupation of electronic states^"—. Both 
mechanisms predict, under appropriate conditions, a neg- 
ative resistivity value in the oscillation minima at a small 
dc current. At some finite current density jo, the resis- 
tivity crosses zero and becomes positive^ i^-'^d^ . The sec- 
ond effect invoked for explanation of ZRS is spontaneous 
formation of current domains with current density jtr^ 
which occurs due to instability of systems with negative 
absolute and/or differential resistivity ^^"^''. In the sim- 
plest geometry the system breaks into two domains car- 
rying equal Hall currents of opposite direction. The Hall 
electric fields related to these currents are also equal in 
magnitude and have opposite directions. The resulting 
dissipative component of the resistivity or conductivity 
tensor is zerc^^. Two works^^'^^ provided an experimen- 
tal evidence supporting the spontaneous domain forma- 
tion in ZRS. 



In this paper we show that a positive diffusion coeffi- 
cient can stabilize a homogeneous state of a finite 2DES 
even though the system has a negative resistivity. We 
point out that combination of a negative resistivity and 
a positive diffusion coefficient naturally arises in 2DES 
under microwave radiation. Solving the transport equa- 
tion jointly with the electrostatic problem for a finite 
2DES, we determine a stability condition which includes 
conductivity, diffusion coefficient, and the system size. 

Further, we allow for a non-zero difference between the 
work functions of two contacts to 2DES (leads), that 
gives rise to the photogalvanic effects in an irradiated 
sampl©i2. We explore the distribution of the electrostatic 
potential and the electron density in such a sample for 
different values of the conductivity. Most remarkably, 
we find that in the diffusion-stabilized regime of negative 
conductivity the system splits into two regions with op- 
posite directions of electric field. We argue that this effect 
is a precursor of the domain structure that is expected 
to emerge spontaneously in microwave-induced ZRS. 



II. MODEL AND BASIC RESULTS 



A. Microscopic conductivity and diffusion 
coefficient 



The local dissipative current density in a weakly non- 
uniform 2DES with surface electron density Ug has a form 



j = -2fjV, 



|z=0 



eDV M^ 



(1) 



where (?!)|z=o is the electrical potential in the 2DES plane 
z = 0, Vs is the surface gradient in this plane, a is the 



dissipative component of the magnetoconductivity tensor 
(per spin), and D is the diffusion coefficient. 

In equihbrium, the conductivity ctq and the diffusion 
coefficient Dq obey the Einstein relation 2cro — e^xDo, 
where x = dug/d^ is the equilibrium static compressibil- 
ity (here /i is the 2DES chemical potential). The Ein- 
stein relation guarantees that j = — 2crVs?7/e is zero in 
the equilibrium state with constant electrochemical po- 
tential 77 = e(?!)|z=o + M; further, the current response does 
not depend on the type of perturbation (electric field or 
density gradient). 

In the presence of microwaves the Einstein relation is 
violatedfii^ 



2cr ^ e\D . 



(2) 



This violation leads to the photogalvanic effects observed 
experimentally in Refs. [l5ll23 and is at the heart of the 
phenomena discussed in the present work. 

We start with demonstration that the theoryii allows 
for coexistence of a negative dissipative conductivity and 
a positive diffusion coefficient. This result is the most 
prominent for the case when effect of the microwave 
radiation on the electron kinetics is governed by mod- 
ification of the electron distribution function (inelastic 
mechanism). ^^ For this mechanism, equations for the 
conductivity a and the diffusion coefficient D read 



a = -cfd 



iy\e) dfje) 
Vn de 



de. 






(3) 



(4) 



Here ajj = ngc'^ /2m*Lu'^T is the Drude conductivity in 
a classically strong magnetic field, uJcT ^ 1; 'ti* is 
the effective mass, Wc the cyclotron frequency, and r 
the transport relaxation time. We remind the reader 
that at uJcT ^ 1 the dissipative resistivity is propor- 
tional to a (since the dominant component of the con- 
ductivity tensor is the Hall conductivity that is only 
weakly affected by microwaves). Further, i/(e) is the den- 
sity of states in disorder-broadened Landau levels and 
i^o — m* /2TTh'^ is the density of states per spin at _B = 0. 
The non-equilibrium distribution function /(e) is deter- 
mined by the kinetic equation>ii Figure [T] illustrates the 
magnetic-field dependence of the dissipative conductivity 
a and of the diffusion coefficient D for typical parame- 
ters. (Details of numerical procedure used for calcula- 
tion of density of states iy{e) and nonequilibrium distri- 
bution function /(e) are given in Ref. l2lh Under mi- 
crowave radiation with the circular frequency w, a shows 
strong magnetooscillations (with periodicity controlled 
by the ratio lu/lJc) and becomes negative around min- 
ima. At the same time, D remains almost unaffected by 
microwaves (apart from Shubnikov-de Haas oscillations 
that are strongly suppressed due to temperature smear- 
ing; for parameters in Fig. [T] their amplitude remains 
within 1%). 
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FIG. 1: Magnetic-field dependence of (a) the dissipative con- 
ductivity CTo, o" and (b) the normalized electron diffusion co- 
efficient D/ Do, shown as a function of the ratio bj^jio, where 
ojc is the electron cyclotron frequency and oj is the circular 
radiation frequency. The quantities ctq and Do correspond 
to the dark conditions, while a and D were calculated for a 
nonequilibrium state under microwave radiation. Parameters 
of calculations: electron surface density ria corresponds to the 
filling factor of the Landau levels A'^ ~ 100 at cJc/i^ = 1, 
kT jYvjj = 0.2, and lot^ — 10. Here T is the dark temperature 
of 2DES and Tq is the quantum scattering time determining 



the Landau level broadening. 



B. Electrostatics and transport in a 2D stripe 

We proceed now by finding a self-consistent solution of 
electrostatic and transport problems in a sample of finite 
length L subjected to microwave radiation and dc elec- 
tric field. The considered geometry is shown in Fig. [5] 
Specifically, we assume that 2DES forms a stripe in the 
(x, y) plane between two infinite metal plates at x = 
and X = L which model Ohmic contacts. The magnetic 
field B is parallel to the z axis. In order to find a station- 
ary spatial distribution of the 2D electron density Ug {x) 
and of the potential 4>{x, z) we use a system of coupled 
equations that includes the transport equation ([1]), the 
continuity equation 



eus 



V. -j^O, 



(5) 
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FIG. 2: Schematic view of the 2DES and metal contacts for 
considered electrostatic problem. 



and the Poisson equation 



eA(f>(x, z) — 47re(7is — ?^+)(5(z). 



(6) 



Here e is the dielectric constant of the medium sur- 
rounding 2DES, and a constant n_|_ represents the frozen 
positive background charge of ionized donors. In the 
bulk of an infinite sample, the electroneutrality requires 
Ug = ?T.+ . In what follows, we assume that relative den- 
sity variations are weak, Ug — n_|_ ^ n_|_, which allows 
us to use spatially- independent a and D in Eq. ^. For 
simplicity, the system is assumed to be infinite in y di- 
rection. This setup can be termed "quasi- Corbino" , with 
X axis corresponding to the radial direction and y to the 
azimuthal direction of the Corbino disk. By symmetry, 
ns{x) and (j){x,z) do not depend on y, and the dissipa- 
tive current ([T]) flows parallel to x-direction. The lo- 
cal density and electric field define the Hall component 
jy{x) = en{x)vu{x), a dissipationless flow in y-direction 
with local drift velocity vd = —{c/B)d(j){x,0)/dx. 

The above equations are supplemented by boundary 
conditions for the electrostatic potential, 

</)(O,z) = 0o, 0(i,^) = O, (7) 

{d(j)/dz),^±oo ^ , (8) 
and for the density, 

ns(L)=n+, ns{0) = n+ - xeHc ■ (9) 

In Eq. ([9]), eUc is the difference of work functions be- 
tween the left lead and 2DES. Due to large density of 
states in the metallic lead, the interface charge density 
in 2DES e{ns{0) — n_|_) is fixed by eUc and is not affected 
by either the radiation or the current fiow (see Ref. [13 for 
details). For simplicity, we assume that the right lead has 
the same work function as 2DES, so that ns{L) = n+. As 
we show below, a finite eUc j^ makes the electrostatic 
problem nontrivial by introducing the spatial variation of 
the density Hg- Integrating Eq. ([1]) over x from to L, we 
get a relation between the current density j, the electro- 
static potential difference (j)o , and the chemical potential 
difference eUc^^ 



jL = 2a(f>o - e\DUc 



(10) 



In the absence of 2DES (ris — n+ = 0) the Poisson equa- 
tion dH]) has a trivial solution (j){x, z) = 4>q{\—x/L). Split- 
ting this contribution and expanding the remaining part 
in eigenfunctions of the Laplace operator with proper 
boundary conditions [Eqs. ([7]) and dH) with 0o =0], we 
get 



(X, z 






Nl/i, 



fe=i 



kirx 



(11) 



With the use of Eq. ©, the charge density can be ex- 
pressed in terms of coefficients Ak as follows: 



e[ng(x) -n+)^ 2_^ —Ak sm — — . 



fc=i 

We integrate now Eq. ([Ij from x to L, 

j{L — x) — 2a(f>{x, 0) -t- eD[ns{x) — n+J 

and substitute here Eqs. pT|) and ([T^ . Using 

X 2 ■r-^ 1 . nkx 
^- Y = -y^k sm — — , 

fe=l 



(12) 



(13) 



(14) 



we finally obtain the coefficients Ak governing the sta- 
tionary solution of the problem: 



A, 



4L e'xUc 



Trek k + L/ttX 



(15) 



Here we employed the relation (|T0)) and introduced the 
non-equilibrium screening length 



A 



eD 

Ana 



(16) 



(discussed in more detail below). 



C. Nonequilibrium screening length 

The nonequilibrium screening length (J16l) which enters 
the solution (|15p . plays a central role in nonequilibrium 
transport in inhomogeneous 2DES. In equilibrium (in the 
absence of microwaves) , the Einstein relation holds and A 
reduces to the conventional 2D Thomas-Fermi screening 
length. 



A^ Ao = 



27re^x' 



2(70 



■xDo- 



(17) 



Using Eqs. ([I}, ([S|), and ^, it is easy to show that 
indeed A replaces the equilibrium Aq in all electrostatic 
problems. Under microwave radiation, the local conduc- 
tivity a in Eq. (|16|) can cross zero and become negative, 
as illustrated in Fig. [T] The screening length A diverges 
at cr = 0. As we show below, in the regime cr < the 
stationary solutions of the linear problem become electri- 
cally unstable on the spatial scale |A| determined by the 
nonequlibrium screening length. 



D. Stability conditions: finite-size effects 

We check now the stabihty of the obtained solution 
with respect to slow spatio-temporal fluctuations. To this 
end, we add a fluctuating part n, (f) and j to n^, cf), and 
j. To satisfy the boundary conditions {n,(l>}\x=o,L — 0, 
we take 



n{x,y,t) — Snq{t)e^'^^'^ sin q^x 



(18) 



with qx — Trk/L, k — 1,2,.., and continuous qy. Using 
Eq. ([6]), the corresponding 

4>{x,y,z,t) - —Sn,{t)e'''yy-'i\'\smqxx, (19) 

where q — {q^ + q^Y^^- The fluctuating part of the cur- 
rent is therefore 



-2aVs4'\z=Q ~ eDVsfi - 2aHeVs 



(20) 



where we added to Eq. ([T]) the Hall term with an 



ensc/ B and txy 
dS]) gives 



^yx 



1. The continuity equation 



— = — A,0 ,=0 + DA,n = -[ + D] q^n. 21 

at e \ eq ' 



(The Hall term drops out since VseVgi^ — 0.) Therefore, 
fluctuations do not grow in time if the stability condition 



A"^+g>0 



(22) 



is satisfied. The condition is most restrictive for soft 
modes with small q. In an infinite system, where pertur- 
bations at arbitrarily long spatial scale q~^ are possible, 
the condition ([^ reduces to the usual one: A > or, 
equivalently, a > 0.^^ In any finite system, q~^ is lim- 
ited by the system size, and instability threshold shifts 
to negative a. In particular, in 2D stripe of width L, the 
minimal wavenumber q — -k / L corresponds to the low- 
est harmonics fc = 1 in Eqs. (fTT|) - (|T5|) . and the stability 
condition ([22) reads 



or, equivalently. 



L/ttA > -1, 



a > -eD/AL. 



(23) 



(24) 



III. ANALYSIS OF RESULTS 

In this section we discuss the obtained result ([15]) for 
the field (ITT|) and density p^ distribution in different 
physical situations. 



A. Homogeneous stable state with negative 
conductivity 

In the plain-capacitor contact configuration. Fig. [21 
and for vanishing difference of the contact work functions, 
Uc ~ 0, all harmonics Ak = 0, see Eq. ([T5)) . The elec- 
tron density p^ in initially homogeneous 2DES remains 
constant, rig = n_|_, independent on external bias eV = 
?7(0) - T]{L) = e(/i(0) - e(/i(L), and 0(x, z) = ¥{1 - x/L), 
see Eq. ([TT|) . The dissipative current ([U reduces to 



J - 2aV/L, 



(25) 



where a manifests the microwave-induced oscillations il- 
lustrated in Fig. [T^. 

In the finite system, the result ([^5)) holds even for neg- 
ative conductivity values, as long as the stability condi- 
tion (|24p is satisfied. Our analysis in Sec. [HD] shows that 
diffusion with D > can stabilize otherwise unstable^ 
homogeneous state with cr < in a finite 2DES. Accord- 
ing to Eq. ([23| . the lowest observable value of negative 
conductivity in a stable homogeneous state is 



-eD/AL ~ -7rAocro/i, 



(26) 



where in the last equality we used D ~ Dq and Eq. (|17l) ; 
the subscript refers to the dark equilibrium state. The 
critical value Oc corresponds to the critical value A^^ = 
—Ti/L of the inverse nonequilibrium screening length (J16p 
fixed by the lowest possible wavevector in the system. 

At cr < (Tc < 0, the system becomes electrically unsta- 
ble and breaks into domains. Possible domain configu- 
rations and their dynamics were addressed in Refs. [23l - 
[26l . The appearance of electric domains implies an ac- 
cumulation of charge at the boundary between domains 
of opposite polarity (domain walls). The nonequilibrium 
screening length sets the spatial scale of modulation of 
electronic density and electrostatic potential in the do- 
main phase. Therefore, a proper account for the violated 
Einstein relation is crucial for understanding the criti- 
cal properties of the nonequilibrium phase transition to 
ZRS, the microscopic structure of the domain walls, their 
dynamics, and sensitivity to boundary conditions and de- 
tails of the disorder potential. 

It should be mentioned that a negative residual con- 
ductivity has also been found in Ref. [2J. Specifically, 
that work considered the domain phase (tr < Cc in our 
terminology) and obtained an exponentially small nega- 
tive residual conductivity. We believe that the exponen- 
tial dependence found in Ref. 23 is most likely an artefact 
of local electrostatic approximation used in Refs. 23 - 261 
Such an approximation is appropriate for description of, 
e.g., the Gunn effect in bulk 3D semiconductors, but is 
not directly applicable to 2DES with quantized motion 
in z-direction. In contrast to 3D geometry, in 2D the 
relation between (j) and Us is non-local, see Eq. ([6]). For 
proper description of the domain phase, one should go 
beyond the linear response in Eq. ([!}. This regime re- 
quires development of adequate approaches for solution 



of arising nonlinear nonlocal equations and will be ad- 
dressed elsewhere. 



B. Photogalvanic effects 

We return now to an inhomogeneous system with a 
difference of work functions in 2DES and the left lead 
eUc < 0. For simplicity, we also assume that temper- 
ature is sufficiently high to suppress the magnetooscil- 
lations of the dark compressibility x ^-nd Shubnikov-de 
Haas oscillations, ^-K^kT jhuJc ^ 1. In this case, D ^ Dq 
(see Fig. [TJd) and x = dug/dfi — ris/fi = ^vq, so that the 
electrochemical potential 77 — e4>\z=o + ns/2vQ. 



1. Current-voltage characteristics 

The knowledge of spatial distributions (|lip and p^ is 
not required for calculation of the current- voltage char- 
acteristics (CVC). Indeed, the CVC is given by the re- 
lation (TTOl) which is obtained directly from Eq. ([T]) us- 
ing boundary conditions ©■ The electrostatic potential 
(^0 = (T^COjO) can be expressed from Eqs. ([7]) and (JH) in 
terms of electrochemical potential drop eV = 77(0) — rj{L) 
(measurable voltage across the sample) a,s (f)o = V +Uc- 
Therefore, the CVC dTOl) reads 



J 



2cr 



V 



2a 



Uc 



1 



A 



(27) 



The last term is responsible for photogalvanic effects. In 
the presence of contact asymmetry Wc 7^ or other source 
of the built-in electric field Uc/L and provided the Ein- 
stein relation is violated, A 7^ Aq, the system displays 
photocurrent j 7^ at zero bias voltage V = Q, and 
photovoltage T^ 7^ at j = 0, see Ref. I17I for details. 
Our present results show that this characteristic stays 
valid for cr < if the condition (|23| is fulfilled. For 
the case of interest, 13 > 0, Eq. ipTl) predicts that both 
the differential resistance dV/dj and the photo-voltage 
V^hoto = Wc(A/Ao — 1) change sign when conductivity 
goes through the zero value. This prediction can be veri- 
fied experimentally, e.g. on narrow Corbino disk samples 
which should behave similarly to the infinite stripe con- 
sidered here. Importantly, the size of the sample L/Xq 
should not be very large in order to enable controllable 
measurements in the stable region ([^ with a < 0. 



2. Field and density distributions 




0.4 0.6 

x/L 



FIG. 3: The electrostatic potential (/>(a;,0)/Wc vs. x/L for 
j — 0, L = SttAo, and several different A/Ao, see Eq. (|31|) . 
Four upper curves correspond to cr > (A > 0). Two lower 
curves are calculated for ct < satisfying the stability condi- 
tion L/ttA > — 1. 



where 



M{X, 



2 v-^ sin kX 2 , -y ^ , s 

- E 1^— = zlm<i>(e'^, 1, >.), (29) 



TT ^ — ' k + X 
fc=l 



and "I>(r, s, u) = J^'iZoi'^ + '^^^ is the Lerch transcen- 
dent. Further, the 3D electrostatic potential reads 



(x,z) 



2LUc ^ sm{TTkx/L) e^'^'^l^l/^ 



TT^An '^^ k + L/ttX 
At z = 0, it can be represented as 



(30) 



,^(x,0) 



V+Uc{l 



A 

a;^ 

+ WcT-AA — ,— 

Ao \ L TTA 



L 



(31) 



Finally, the electrochemical potential 



77(x) - r,{L) =eV[l-- 



^eMc 1 



-N 



TTX L 



(32) 



The charge and field distribution in the interior of the 
sample provides an additional insight into the problem 
and uncovers an interesting behavior in the vicinity of 
the instability threshold. Using Eqs. ([n|)-([15l) and (P7)) 
we obtain the density profile 

n,{x) - n+ _ fi{x) - fi{L) _ ,^ , -^ ^ , ,^^. 



ns{0) 



eUc 



,., TTX L 



3. Enhancement of built-in field 

Assuming an open circuit, j — 0, in Fig. [3] we show 
how the built-in electrostatic potential (j){x, 0) is modified 
through the microwave-induced changes of the screening 
length A. Equation (P7|) for the case j — yields 



0(0, 0)/Uc = V/U, + 1 = A/Ao ~ ao/a . 



(33) 



In the last equality, we used D ~ Dq. The electro- 
static potential drop along the sample is enhanced at 
< (T < (To (i-e. A/Ao > 1), which is illustrated by two 
upper (green) curves in Fig. [3] (the colors in Fig. [3] and 
|4] refer to the online version). The black curve (third 
from above) illustrates the equilibrium distribution. In 
the case u > ctq (orange curve A/Aq — 0.5 in Fig. [3|), the 
potential drop diminishes with respect to the equilibrium 
case. This behavior at cr > illustrates the enhancement 
(cr < (To) and suppression (a > ao) of the built-in electric 
field proposed for interpretation of photogalvanic effects 
in Ref. |22 on the experimental basis. 

At (T — >■ 0, the potential cJ)q = (/)(0, 0) diverges in 
an open circuit, (namely, (j^o ~^ +oo at cr — )■ -1-0 and 
(j)Q — >■ — oo at (T — ;■ —0). This divergence is cut by nonlin- 
ear corrections to a and D which are not taken into ac- 
count here.^^ At cr < and j = 0, illustrated by two low- 
est (red) curves in Fig.[3l the potential (j>{x, 0) develops a 
single minimum which approaches the position x ~ L/2 
when a approaches the critical value <Tc = —TrXoao/L, 
Eq. ([26t . In Fig.|3]we took L = SttAo, so the correspond- 
ing critical value of A is Ac = —L/n = — 3Ao. 



4-. Analysis of the field and density distributions in 
different regimes 

In an open circuit, j — 0, illustrated in Fig. [3l lin- 
ear terms oc (1 — x/L) in ()31|) and (|32l) cancel out [see 
Eq. (HZl)]. Therefore, 



fi{x) - fijL) 
-eUc 

^ -nix) - TiiL) 

(A/Ao - l)eUc 






^^(^'«) 



i = o. 



(34) ^ 



General behavior (including j ^ 0) is illustrated in Fig. SI 
where we subtract the linear term oc V{1 — x/ L) in 
Eqs. (PT|) and (|5^ . The rest, expressed in units Wc, is 
fully determined by two parameters L/ttA and A/Ao. The 
relation between j and V for an arbitrary measurement 
scheme follows from CVC (P7)) . 

Several parametric regions deserve special attention. 

(i) Continuous limit \ <^ L. This limit of strong 
screening is illustrated by three upper curves (green and 
black) in Figs. 0^, b, and c. In the region a;, A <C i, the 
summation in Eq. (j29p can be replaced by integration. 
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OO 

2 r dt 



TT J t 

x/X 

2 



— ci— sin si— cos — 

TT A A TT A A 



A,x<L. (35) 



Here si and ci are the integral sine and cosine defined by 
ci(z) -|-i si(z) = — /. (e**/^) dt. Similar distribution was 
obtained in Ref. |28| for equilibrium semi-infinite 2DES 
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FIG. 4: Spatial distribution of (a) charge density, [/i(a;) — 
^{L)\/e'Uc vs. x/L, Eq. ((28]); (b) the electrostatic potential, 
[<t){x,Q) - V{1 ~ x/L)\/Uc vs. x/L, Eq. dMl; and (c) the 
electrochemical potential, [rj{x) — r](L) — eV{l — x/ L)]/eUc 
vs. x/L, Eq. ^. In (b) and (c) L = SttAq as in Fig. El The 
parameter L/ttA = {100, 10, 3, 1, 0, -0.5, -0.7, -0.8} from top 
to bottom. 



with different boundary conditions at i? = 0. At large 
distances, Eq. (j35]) produces a power-law decay. 



N ~ 2A/7rx , 



A<a; <L, 



(36) 



characteristic for 2D electrostatics. 

(ii) Enhanced screening, A ^ Ao, see two upper (green) 
curves in Figs. |4^, b, and c. This case corresponds to the 
well developed maxima of the microwave-induced resis- 
tance oscillations, a ^ a^. Beyond the screening length. 



X > X, both and t] change linearly with x. The slope, 



V,77 



= -V. 



(x,0) = ^ = -, (37) 



is given by Eqs. ^^ and 1^ with A/" -J> 0. The second 
(diffusion) term in Eq. ([T]) is relevant only in the nar- 
row region x < A, where |Vs/i/e| ~ \'Vs<p\z=o\Xo/X ^ 

(iii) Equilibrium, A = Aq, see black curves (third from 
top) in Figs. H^, b, and c. When Uc ^ but there is no 
microwave field (A — Aq), (j){x,0) and ns{x) vary near 
the left contact on scale Aq, while the nonlinear part 
of r]{x) is zero. Since the Einstein relation holds, two 
terms in Eq. ([T|) combine into one, j — — 2crVs77/e = 
~2ei>QD\7 sTj — 2aoV/L. In other words, two terms 
in Eq. ([ij partially compensate each other such that 
V^Ty = —eV/L remains constant in space despite both 
Vs?is and Vget/) vary strongly with a; < Aq. 

(iv) Diffusion dominated screening, A > Aq, see or- 
ange curves (fourth from top) in Figs. |4^, b, and c. In 
this regime, diffusion dominates in the sense e^xD > 2a. 
The nonequilibrium screening charge distribution [which 
is smooth compared to the equilibrium case (iii)] pro- 
duces an "overshoot" in (/)(x, 0): two regions with the op- 
posite orientation of induced electric field appear, Fig.|3lD 
(the total field may not change sign if the external voltage 
is sufficiently large). The corresponding nonequilibrium 
correction to 77, Fig. ^t, changes sign and remains nega- 
tive at A^^ < Aq"^ (including negative A~^, see below). 

(v) Zero-conductivity state, A~^ = 0, see blue curves 
(fourth from bottom) in Figs. |4^, b, and c. In the limit 
L — > 00, homogeneous state of 2DES becomes electri- 
cally unstable at A~^ = OJ^ In the finite system, the 
instability threshold shifts to negative a determined by 
the condition L/ttA = — 1. 

As was discussed in Sec. IIII B 3[ in an open circuit, 
j = 0, the electrostatic potential 0(0, 0) has a singu- 
larity at CT = 0, which implies the necessity to include 
the non-linear effects. This divergence does not appear 
if one fixes the voltage V instead of current. Indeed, 
Eq. dD) with cr = and boundary condition © yields 
the charge density varying linearly with a;. A/" = 1 — x/ L, 
see Fig. |4^. The corresponding current density ([T|) has a 
F-independent value 



js = eD{no - n+)/L = -xe'^DUc 



(38) 



fixed by the boundary conditions. The voltage V de- 
couples and can be arbitrary within the range where the 
linear-response approximation ([1]) is justified. 

At A~^ = 0, second term in Eq. ([5^ (which contains an 
indeterminate form of the type cx) • 0) can be represented 



77(x)-r?(L)-ey(l-f) _ 2L 



eUc 



TT^XO 



ImLi2(e'""/^), (39) 



see Fig. gj:. Here we used X^fcli ^ ^ ^i^i T^ = 
Im Li2(exp[i7ra;/L]), where Li2(z) is the dilogarithm func- 



tion. The potential profile induced by the linear vari- 
ation of charge density A/" = 1 — x/L is e(j){x,0) = 
r]ix) - 7]{L) + eUcil - x/L), Fig. Hb. 

(vi) Stable negative conductivity state, —n/L < X~^ < 
0, see three lowest (red) curves in Figs. [3^, b, and c. The 
distributions ([^S]) . ([5T|) . and ([5^ are dominated by first 
harmonics in Eq. (P^j) . 



J\f 



TTX L 



2 sin(7rx/L) 



(40) 



L ' ttXJ tt I + L/ttA' 
which diverges at the instability threshold, l + L/vrA = 0. 



5. Towards the domain structure 

The divergence of our solution at the threshold A = 
— L/tt signals an instability and transition to the domain 
phase. We believe that inclusion of nonlinear effects (in 
particular, taking into account dependence of a entering 
Eq. ([1]) on the electric field V s4>) should make the theory 
applicable also in the domain regime 1 + L/tiX < 0. Work 
in this direction is currently underway. 

Let us emphasize the emergence of a nonmonotonous 
profile of the electrostatic potential [see Figs. 3 and 4(b)] 
implying formation of two regions with opposite direc- 
tions of the electric field. In our solution, the direction of 
electric field (H — vs. — h) in these two regions is deter- 
mined by the sign oiUc- A non-zero U^ explicitly breaks 
the inversion symmetry L/2 + x ^ L/2 — x. We argue 
that this effect which emerges at A < and grows while 
the system approaches the instability threshold, is the 
precursor of the domain structure which fully develops 
at A^^ < —tt/L. As was discussed in Sec. IIII A| in the 
homogeneous case Uc = the inversion symmetry is pre- 
served at (Tc < cr < and gets spontaneously broken at 
the instability threshold CTc — o. 

A further important direction for future work is to go 
beyond the mean field approximation and to study the 
effects of fluctuations and noise on the transition into 
the ZRS regime. It is expected that field-theoretical ap- 
proaches developed for nonequilibrium phase transitions, 
in particular, Martin-Siggia-Rose formalism, will be use- 
ful in this respect. First steps in this direction were made 
in Ref. [H 



IV. SUMMARY 

In summary, the effect of the microwave radiation 
on the electron energy distribution function of a 2DES 
causes giant magneto-oscillations of the conductivity rel- 
ative to its dark value and practically does not alter 
the electron diffusion coefficient. Such effect leads to 
magneto-oscillations of the photo-galvanic signals and of 
the screening length, which affects the potential profile 
in a sample. At the oscillation minima, the conductivity 
can become negative, which leads to a negative value of 
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the non-equilibrium screening length. We have derived 
the stability condition at which a finite 2DES can pos- 
sess a stable state with a negative conductivity. When 
the conductivity becomes negative, the differential re- 
sistance and the photo-voltage also change their signs, 
which can be observed experimentally. We have further 
solved the combined transport and electrostatic problem 
and determined the profiles of the potential and the elec- 
tron density inside the sample. In the stable state with 
a negative conductivity, the potential profile consists of 



two regions with opposite directions of the electric field. 
The amplitude of these fields increases when the system 
approaches the instability threshold. This effect is a pre- 
cursor of the domain structure in the regime of sponta- 
neously broken symmetry. 
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